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Abstract. Radial sonic crystals (RSC) are fluidlike structures infinitely periodic 
along the radial direction. They have been recently introduced and are only possible 
thanks to the anisotropy of specially designed acoustic metamaterials [see Phys. 
Rev. Lett. 103 064301 (2009)]. We present here a comprehensive analysis of two- 
dimensional RSC shells, vifhich consist of a cavity defect centered at the origin of 
the crystal and a finite thickness crystal shell surrounded by a fluidlike background. 
We develop analytic expressions demonstrating that, like for other type of crystals 
(photonic or phononic) with defects, these shells contain Fabry-Perot like resonances 
and strongly localized modes. The results are completely general and can be extended 
to three dimensional acoustic structures and to their photonic counterparts, the radial 
photonic crystals. 
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1. Introduction 

Acoustic metamaterials are a new type of structures with exciting properties. They 
consist of periodic arrangements of sonic scatterers embedded in a fluid or a gas and their 
extraordinary properties appear at wavelengths much larger than the lattice separation. 
Very recently, these authors employed acoustic metamaterials to demonstrate that is 
possible to create a new type of sonic crystals named Radial Sonic Crystals (RSC)[1]. 
These crystals are radially periodic structures that, like the crystals with periodicity 
along the Cartesian axis, verified the Bloch theorem and consequently, the sound 
propagation is only allowed for certain frequency bands. Some interesting applications 
of RSC have been envisaged. For example, they are proposed to fabricate dynamically 
orientable antennas or acoustic devices for beam forming[lJ. 

In this paper we present a comprehensive analysis of RSC shells, which are 
structures consisting of a central cavity surrounded by finite RSC, both being embedded 
in a fluid or a gas. It will be shown that these structures generate two types of resonant 
modes; ones are strongly localized in the cavity the the others are weakly localized in 
the surrounding shell. Particularly, those strongly localized or cavity modes can have 
extremely large quality factors and made them suitable for developing acoustic devices 
in which localization is specially needed; for example, it can be used to enhance the 
phonon-photon interaction in mixed phononic-photonic structures. Moreover, they can 
be employed as acoustic resonators of high Q (quality factor) since it is shown here 
that acoustic cavities based on RSC shells can achieve Q-values much larger than that 
reported up to date[2]. 

The paper is organized as follows. After this introduction, in Section 2, the theory 
behind the infinitely periodic RSC is briefly introduced for the comprehensiveness of 
the work. Section 3 studies the resonance properties of a particular RSC shell and 
special emphasis is put in getting analytical formulas that give a physical insight of the 
properties associated to the two type of resonant modes. The mode pressure patterns 
are analyzed in Section 4. Finally, the work is summarized in Section |5l where we also 
give a perspective of possible application of these structures as well as an account of 
future work to be performed. 

2. Radial Sonic Crystals 

The propagation of a harmonic acoustic wave, P{r;u), in an anisotropic fluid is given 
by the following wave equation [3\: 



where B and p ^ are the bulk modulus and the reciprocal of the mass density tensor, 
respectively. 




(1) 
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By assuming that the propagation takes place in 2D, the pressure field P can be 
factorized as follows, 

P{r,9;u) = J2P,{r)e^''\ (2) 
where {r,6) are the cylindrical coordinates and Pq{r) accomplishes the equation 

P,ir) = 0, (3) 
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where q are discrete number (g =0,1,2, ...)pr and pg being, respectively, the components 
of the anisotropic mass density tensor in 2D. 

Due to anisotropy, all the coefficients in ([3]) are independent and, therefore, equation 
([3]) can be made invariant under the set of translations r — )■ r + nd, where n is an 
integer and d is the period, which is prefixed parameter. Acoustic structures having 
these properties are called Radial Sonic Crystal (RSC)[1] because they are formally 
equivalent to the crystals already known, like atomic crystals, photonic crystals or 
phononic crystals, whose corresponding wave equations are invariant under translations 
along the Cartesian axis. 

2.1. Acoustic parameters radially periodic 

The conditions that leave ([3]) invariant under translations of the type r r + nd are: 

r + nd r , , , 

(4a) 



Pr{r + nd) p,.{r) 
r + nd r 



B{r + nd) B{r) 
{r + nd) pe{r + nd) = r po{r) (4c) 

From these conditions is easily concluded that the acoustic parameters can be cast into 
the formfl]: 



pr{r) =rpr{r) (5 a) 

Pe\r)=rpg\r) {5b) 
B{r) =rB{r), (5c) 

where Pr{r),pg{r) and B{r) represent periodic functions of r. Let us remark that any 
periodic function is valid to get the invariance obtained from (4). 

Figured] depict, as a typical example, the profile of Pr{r), which has been generated 
by the product of r with a radially periodic function p^, as explained in (5a). Notice 
that the RSC is only defined in the region 2< r/d < 10. So, this plot also illustrates 
the case of an acoustic cavity generated by a RSC shell of 8 layers; below r/d =2 the 
mass density is pa and for r/d >10 the surrounded background has density pb- 
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Figure 1. Plot of a radially periodic parameter. It is seen how there is a divergence 
at the origin, but in practice the finite structure will consist in a cavity region, the 
crystal and the background. 



2.2. Wave Propagation and Band Structure 

The periodicity of coefficients (12. ip allows to apply the Bloch's theorem to solve (|3]). 
Thus, the solutions Pq{r) can be chosen to have the form of a plane wave times a function 
with the periodicity of the lattice[l]: 

nW = e'''''E V"' (6) 



n 



where kq is the wavenumber and Gn are the reciprocal lattice vectors given by Gn = 
2n7r/d, n being integers. 

The periodicity of coefficients allows the expansion in Fourier series: 

Pr n 

— = T.9ne'''''' (76) 
rpe n 

^ = E /^ne^^- (7c) 

The insertion of these expressions in ([3]) and after using the expansion of the field (E]), 
the following matrix equation is found 

MP = u^NP, (8) 

where the matrix elements are 

Mmn = fm-n{Kg + G^){Kg + G^) + 9 m—n 

(9 a) 
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Figure 2. Acoustic band structure for a two dimensional radial sonic crystal. 

This matrix equation is a generalized eigenvalue problem that can be solved 
numerically with a wide variety of methods. Thus, once the acoustic parameters 
defining the crystal are given, we can solve the eigenvalue equation u = u^Kg), 
obtaining a dispersion relation that is different for each multipole q and noting that 
u{K,) = u{K.,m. 

The final form for the field inside the crystal is given by a linear combination of 
plane waves. 

P(r, e) = Y, e^^^^e^"^ ^<?ne'^" (10) 

q n 

which is a linear combination of plane waves. 

In [I], we have considered a RSC made of two types of alternating "homogeneous" 
layers. Actually, the layer are not exactly "homogeneous" because they need to be 
radially dependent. The term "homogeneous" applies here because this medium is the 
equivalent to the system "ABABABA..." employed in standard heterostructures. The 
difference is that in this case Pr,Pd and B are considered as the "stepped" functions. 

2.3. Low Frequency Bandgaps 

Figure [2] shows that in the low frequency limit only the dispersion relation corresponding 
to modes q =0 (black lines), has a linear behavior near the zero frequency. The bands 
corresponding to higher modes, q >0 present low frequency bandgaps. In other words, 
for non monopolar modes, the RSC behaves like a "radially homogeneous" medium; 
that is, a system where the periodic functions pr, pe and B take constant values p^, JJg 
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and B, respectively. Then, in such a medium the wave equation Q becomes 

P,(r) = (11) 



B ^ Wo 



The solutions of this equation are plane waves 

P,(r) = e''^\ (12) 

where 

kl = u%-q'~P^, g = 0,l,2,... (13) 
^ B pe 

Therefore, below certain cut-off frequencies given by 

B 



\ 



(14) 



the wavenumbers take imaginary values and evanescent waves are the only possible 
solutions to ([3]). The cut-off frequency defines the low- frequency bandgaps of 
the corresponding mode and it is a property characterizing the so called "radially 
homogeneous" medium, whose acoustic parameters (density and bulk modulus) linearly 
increases with the radial distance to the origin of the crystal. 

3. Resonances in a RSC shell 

Like with other type of crystals, it is not possible to handle infinite RSC. So, let us 
consider a RSC shell consisting of 8 periods as the one described in Fig. [H The shell 
has an inner cavity of radius Ra = 2d that acts like a defect. A scheme of the structure 
under study is shown in figure [H, where it can be seen that the cavity medium A is 
enclosed by the anisotropic shell S and both being surrounded by a background fluid 
B that exists for r > Rf,. Both the medium A and B are homogeneous and isotropic 
fluidlike materials with acoustic parameters pa, Ba and p^, Bh, respectively. The shell 
is assumed to be an anisotropic fluid-like material satisfying |3l 

Here, we develop a formalism similar to the transfer matrix method ^ to study two 
type of scattering problems. First, we consider the scattering of an external sound wave 
impinging the RSC shell. Afterwards, we analyze the case when the exciting sound is 
inside the cavity shell. 

In both media, A and B, the fields are given by a linear combination of Hankel and 
Bessel functions, 

P(r, e) = Y: [Ct,H,{k,r) + Ci^Uhr)] e^^^ (15) 
q 

where ki = uj/ci and i = a,b. 

For the anisotropic shell S, which covers the region Ra ^ r < R^, the pressure field 
can be expressed as 

P(r, e) = Y: [C+0+(r) + C-0-(r)] e^^^ (16) 
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where 0+ and 0+ are two hnearly independent solutions of the wave equation. The 
proposed solutions in the shell are completely general and it is assumed that the acoustic 
parameters pr{r), pg{r) and B{r) are arbitrary functions of the radial coordinate r. Later 
we will consider a particular RSC. 

To determine the relationships between coefficients for £ = a,h,s boundary 
conditions must be applied at the boundaries Ra and R},. The boundary conditions are: 
(l)the continuity of pressure and (2) the continuity of the radial component of particle 
velocity. Because of the radial symmetry of the problem it is possible to uncouple the 
equations for the different multipole components q. 

Therefore, at r = i?^, for ^ = a, 6, we require that: 

C+H,{kM + Ci^UhR,) = Ct^MRt) + C-0-(^€) (17a) 



C/MkeRe) + CiXikRi)] = \CtAMr) + C," (r)] ^^^^ (176) 



Pi ^ Pri-n-^J 
Note that, for convenience, the continuity equation of particle velocity has been multiply 
by the radial coordinate Ri at both sides. 

Solving for the coefficients Cf^ the transfer matrix M of the shell is found 

After some algebra, the matrix M can be factorized as M = M~^M5Mb, where, 

^aRaPa JqiJ^aRa) Jqi^aRa) \ {AQ \ 

« ~^[-kaRap-'H'^{kaRa) H,{kaRa) J ^ ' 

M = ( '^'^ Rb^PriRb)(f)l2 \ ..g,^ 

V RaP;\Ra)hl R^'RaPr{Rb)p;\Ra)h2 J ^ ' 

^ f H,{kM j,{kM \ 

' \ hRbPb'H'^ihRb) kbRbPb'J'gihRb) J ^ ' 

The quantities (pij in Mg are Wronskians of the shell function (jy^ir) and their 
derivatives, so that no analytical simplification can be done here. However it is easy tho 
show that the determinant of is one (sec the Appendix). 

Now, the cocfiicients for the two type of scattering problems of interest here can be 
obtained by applying the transfer matrix formalism. 

First, let us consider the case when there is no sound source inside the the cavity 
and, therefore, C+ = 0. The exciting field is an external wave defined by = and 
the goal is to find the coefficients = T'^A^^ and = T^'^A^^, by solving 

/ \ / Mil 



M21 M; 




(20) 



That is. 



Ml 9 

T''^^ (21a) 



Mil 

Mil 



rj^in ^ M22Mii-Mi2M2l ^ ^215) 
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where the numerator of T^" is the determinant of M, which is |M| = |Ma | IM^I |M;,| = 
Pa/Pb- This yields 

Ml 9 

Tf = -'-^ (22a) 
"Mil 

™ ^P^J_ f22b) 
" PbMii ^ ' 

If the exciting field is put inside the cavity, the pressure wave is defined by the 

coefficients C^^ = A^. Now, a standing wave is excited inside the cavity = RqA'^ 

and an acoustic field defined by = TgA^ and = is radiated by the cavity. 

This problem is solved by obtaining the coefficients Rg and Tg in the following matrix 

equation: 




(23) 



Again, solving for the coefficients, 

= ik 

For the two scattering problems analyzed above, the behavior of their corresponding 
coefficients is controlled by a common denominator; the diagonal element Mu. 
Moreover, T^*" = Tg when both the cavity and the background are the same medium; 
i.e.. Pa = Pb- 

If the shell is made of a RSC with A^-periods, the matrix of each period is 
identical, due to the periodicity defining the RSC. The transfer matrix of such system 
is therefore 

M = m-^m^mb (25) 

Since the matrix has modulus one, its eigenvalues are of the form A± = e^*^?'^, 
being Kg = Kg{u) the dispersion relation that defines the band structure of the crystal. 
These properties make possible express the A^-th power of the matrix Mg as [5] 

^^sinNKgd ^sm(N-l)Kgd , , 

smilgd smAqd 
and the transfer matrix of the system as 

1 sinA^A:„d , ^ , ^ sin(A^ - l)A:„ci 

sm Kgd sm Kgd 

The complex frequencies u solving equation Mu = completely characterize the 
resonant modes in the structure; i.e., their frequencies and lifetimes. However, the 
numerical solution of such equation is a cumbersome task and it is beyond the scope 
of the present work. Here, we prefer to give a physical insight of the resonances u for 
which the coefficients have large or small values. 

Two types of resonances can be distinguished: cavity modes and Fabry-Perot like 
resonances. The former are strongly localized inside the cavity shell, the later are due 
to the finite thickness of the RSC making the shell. 
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Figure 3. Localized modes and Fabry-Perot resonances for a finite RSC consisting of 
8 periods and an inner cavity of radius Ra = 2d (see Fig. 1). Top panels: Graphic 
resolution of Eq. (20) for modes and 1, respectively. The dots define the frequency 
positions at which the mode cavities are predicted. The shadow zone define the acoustic 
bandgaps shown in Fig. 2. Bottom panels: The corresponding scattering coefficients, 
To and Ti, which have peaks at the frequencies approximately determined by Eq. (20) 
(see text). 



3.1. Cavity modes 

These resonances occur at frequencies within the bandgaps of the RSC. They are 
characterized by their strong locahzation inside the cavity shell as it is shown below. 

Bloch wavenumbers within the bandgaps are complex number, Kg ^ iKg, and as 
a consequence the periodic functions of NKqd in (j3]) become hyperbolic functions of 
NKqd. Therefore, the dominant term in the expression for Mn are proportional to 
sin NKqd ~ ^e^^'^'^. If we consider, for example, that the sound source is placed inside 
the cavity shell, the field radiated out of the shell is given by (24b), that is 

Tq = 1/Mn ^ e-^^'^ (28) 

This exponentially decreasing behavior with increasing N has been previously 
characterized in other type of crystals, like photonic or phononic crystals. For example, 
the attenuation by a phononic crystal slab has been demonstrate to be exponentially 
decreasing with the number of crystal periods|6]. 

However it is also possible to demonstrate that, in the bandgap region, Tq have 
strong transmission peaks, which are associated to the existence of cavity modes. Thus, 
after some algebra, the Mu element can be cast as 



11 



XqHqihRb) + YqhRbP^'H'ihRb) (29) 
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Figure 4. The same representation than in Fig. |3] but for modes q —2 and q =3, 
respectively. Notice that the large frequency bandgaps (shadow regions) appearing 
at low frequencies (see also Fig. [2] The Fabry-Perot like resonances are shown in the 
narrow passband frequencies (white regions) . 



where 

X, = M^,Kp-'RaJ'g{kaRa) ' UKRa)M^, (30a) 
and 

Y, = M^^kap-'RaJ',{kaRa) ' J,{kaRa)Mg (306) 

When = we have that 

Y, = -UkaRa)/M^, ^ e-^^^^ (31) 

which means that now Mn behaves inversely as before; i.e., the transmission coefficient 
increases exponentially with the number of layers. 

Tg = l/Mii ^ e^^^'^ (32) 

A similar behavior is found after analyzing the case = 0. 

The resonance frequency can be obtained from the condition = 0. The following 
condition must be hold 

, Ip J'qikqRa) _ M^^ _ , ^21 ^ooN 

^"^-^ ^IjKRa) ~Wl^ ^ ^ ^ 

where the quantities (pij have been already introduced as factors in the elements of 
Ms [see (19b)]. Since we are working at bangap frequencies, the ratio 02i/0ii can be 
approximated by 

'^^iK, = -Im[K,{uj)] (34) 
011 
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Figure 5. (a) Coefficients Ti of a RSC sliell witli tlie same parameters tlian in Fig. 
[3l Nut now the calculation is performed with frequency steps of 10^'' units, (b) 
Coefficients T2 of a RSC shell with the same parameters than in Fig. H) But now the 
calculation has been performed with frequency steps of 10^^^ units. The shadowed 
regions in both figures define the low frequency bandgaps of the dispersion relations 
with symmetry q —1 and q =2, respectively, in the corresponding RSC (see Fig. [2]). 



Therefore, the equation for the frequency ujq of the cavity modes is rapidly find 

LOoJ'n{uJoRa/Ca) = ^^^^Im[Kg{uJo)]Jq{UoRa/Ca), (35) 

which can be solved graphically. 

Figures [3] and m show the graphical solution of fl35|) . where the left hand side is 
depicted with black dots symbols while the red lines represent the right hand side. The 
color dots at the crossing points give the (approximate) values of the cavity modes. 
These values are confirmed by the numerical solution of coefficients Rq and T^, which 
are depicted in the bottom panels of Figs. [3] and HI Note how, at frequencies within the 
bandgaps (shadowed zones in the top panels), the coefficients show peaks that practically 
coincide with the points giving the solutions of ( l35l) . The peaks in the passband regions 
(white zones) correspond to Fabry-Perot like modes that are discussed in the following 
subsection. Let us remark that coefficients Rq and Tq in Figs. [3] and H] have been 
calculated with a resolution of 10~^ along the frequency axis and, therefore, their values 
at the peaks are not necessarily converged. For an exact determination of peak heights 
lower frequency steps should be chosen till the convergence is obtained. 

The behavior of coefficients Tq near resonances suggests that their moduli have the 




Figure 6. The quality factor of three different modes localized in the cavity of a RSC 
shell made of N periods. The mode with frequency 0.1973 has dipolar symmetry {q =1) 
while modes with frequencies 0.3167 and 0.5538 both have quadrupolar symmetry 

(q =2). 



following functional form 

\TJuj)\^ = 1 (36) 

The quality factor Q = uq/Au, where Au is the full width at half maximum (FWWHM) 
of the resonance profile. In other words, Au = w+ — where u± being the frequencies 
at which |Tg(co'-|-)p = |Tg(wo)P/2. It is easy to show that 

IB 

uj± = uJo±J — , (37) 

and therefore 




g = 2W^e^^''^ (38) 



This formula predicts that the Q—ia.ctoT of a cavity mode grows up exponentially 
with the number of layers that form the RSC shell. This result is also known in 
the field of photonic crystals, where it was numerically demonstrated, for example, that 
the Q— factor of a cylindrical defect surrounded by a finite metallic 2D photonic crystal 
exponentially increases with the number of layers surrounding the cavity defect [7]. 

To demonstrate the validity of fl38|) we have calculated Q as a function of N for 
three different cavity modes. It must be pointed out that Q is obtained directly from the 
peaks in Tq and, therefore, a very accurate determination of peak's height at the resonant 
frequencies is needed. This requirement is accomplished by exploring the frequency axis 
in very small steps until convergence is confirmed. As an example. Fig. [5] shows the 
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calculation of Ti and T2, which were already depicted in [3] and SJ by using frequency 
steps of 10"'^ and 10~^^, respectively. Now, Tg have values at the frequency resonances 
much larger than in Figs. |3] and H] since are converged. 

Figure [6] represents the Q— factor calculated for several values of N. We have 
analyzed one dipolar mode and two quadrupolar modes. The dipolar mode, which 
has frequency 0.1973, corresponds to the resonant peak having the lowest frequency 
appearing in Fig. |5t^a). The two quadrupolar modes have frequencies 0.3167 and 
0.5538, respectively, and are shown in Figj5]^b). It is observed that the predicted 
exponential dependence is confirmed by the numerical simulations and provide the order 
of magnitude of Q for these new structures. These structures used as acoustic resonator 
could provide Q values much larger than that already obtained by perfect cylinders or 
spherical resonators |2] 

3.2. Fabry-Perot Resonances 

Fabry- Perot like resonances are generated by the oscillating terms sinNKgd and 
sin(iV — l)Kqd in where Nd defines the RSC thickness, Nd = — Ra. These 
resonances are confined inside the RSC region. They are equivalent to Fabry-Perot 
modes already characterized for sonic crystals slabs [8j. A typical feature of these 
resonances is that they are frequency equidistant in the regions where the dispersion 
relation is linear. This feature is clearly observed in the profile of coefficients R^ and 
To shown in Fig. [3l where the peaks in the lower part of the first frequency band are 
equally spaced. 

The Fabry-Perot frequency oscillations, which are produced by the periodic 
functions of Kg{Rb — Ra) in (l27|l . are coupled with the oscillating behavior of the quasi- 
periodic functions with argument kaRa embedded in M^. Therefore, two parametric 
regimes can be distinguished. First, the regime of thick shells occur when Rh » Ra 
and, consequently, resonances localized in the RSC are dominants and the cavity modes 
will be scarcely shown and will appear as additional peaks on top of the oscillating profile 
of Rg and Tg. The second regime is define by the condition Rf, ^ Ra and correspond to 
the case of thin shells. In this case the oscillating profile will be dominate by resonant 
peaks associated to cavity modes and, therefore, Fabry-Perot oscillations, if any, will 
appear in between. 

Fabry-Perot modes have Q— factors much lower than that of cavity modes. For 
example, for the mode with dipolar symmetry located at frequency 0.3649 [see the arrow 
in Fig. Et^a)] we get Q =380. The modes with quadrupolar symmetry located at 0.5936 
and 0.6122 [see the arrows in Fig. Mjo)] have Q values of 4151 and 3364, respectively. 
Lower Q— values means less localization and, consequently, easy excitation by external 
sources. At this point it is interesting to point out that the existence of a great number 
of resonances combined with their easy excitation could be employed to design acoustic 
devices for broadband detection of external sound sources as was suggested in [Ij. 
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Figure 7. Pressure maps (amplitude) of the cavity mode with frequency 0.1973 (in 
reduced units) and dipolar symmetry existing in the RSC shell made of 10 layers (see 
also Fig. |3]). (a) Mode excited by a punctual sound source (not depicted) put at 
position {x,y) =(0,1) inside the cavity, (a) Mode excited by an external sound wave 
with a plane wavefront. The white circles define the borders of the RSC shell. 
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Figure 8. Pressure maps (amplitude) of the two cavity modes existing in the RSC 
shell with quadrupole symmetry (q =2). Their frequencies are 0.3159 (a) y 0.5564 (b) 
in reduced units. They are excited by a punctual sound source (not depicted) put at 
position {x,y) =(0,1) inside the cavity. The white circles define the borders of the 
RSC shell. 
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Figure 9. Pressure maps (amplitude) calculated for the scattering of a sound plane 
wave impinging a RSC shell made of 8 layers, (a) External sound with frequency 
0.3159 (in reduced units), which is resonant with the first quadrupole cavity mode 
in Fig. EKb). (b) External sound with frequency 0.5564 (in reduced units), which is 
resonant with the second quadrupole cavity mode in Fig. [5{b). 



4. Modes pressure fields 

In this section we use pressure maps to discuss in brief tlie salient features of the different 
resonant modes that can be excited in a RSC shell. As typical examples we consider 
modes with dipole (g =1) and quadrupole (g =2) symmetry to illustrate their behavior 
under different excitation conditions. Cavity modes and Fabry-Perot like modes are 
shown. As it is pointed out before, the main feature distinguishing the cavity modes 
from the Fabry-Perot like modes is their difficulty of being excited by using external 
sound sources. 

Figures [Tl^a) and UIJo) depict a cavity mode with dipole symmetry localized in a 
RSC shell made of 10 layers. The maps are obtained under the two different excitation 
conditions reported in Section 3: (a) by using a punctual sound source placed out of 
the cavity center and (b) by using an external sound with a plane wavefront. In both 
cases the exciting sound has the same frequency than the excited mode. Note the in 
both cases the cavity mode is excited and strong localization inside the cavity is clearly 
observed. The black circle in [Tl^a) represent the region not calculated by our computer 
code due to the presence of the exciting source. 

Figures M and M represent the pressure maps of the two quadrupolar modes whose 
Q-factors were analyzed in Fig. |6l Figure [8] has been obtained by using an inner 
punctual source while Fig. |9]was obtained by using the scattering of an external sound 
wave interacting resonantly with the cavity mode. In [9] the thickness of the RSC shell 
is 8 layers since the modes are easier excited from outside. The oscillations observed in 
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Figure 10. Pressure maps (amplitude) of the Fabry-Perot mode of dipole symmetry 
{q —1) that appears at frequency 0.3648 (see Fig. [5] (a) Excitation by a punctual sound 
source (not depicted) put inside the central cavity; at {x,y) = {0,d). (b) Excitation 
by external sound with a plane wavefront. The black circles define the borders of the 
RSC sheh. 



m^b) out of the shell are produced by the external sound coming from the left side of 
the X— axis. Like the cavity mode with dipole symmetry, the quadrupole modes are also 
strongly localized inside the cavity shell. 

Finally, Fabry-Perot like mode with dipole and quadrupole symmetry are shown in 
Figs. [TOl [11] and [12] under the two different excitation conditions. Note that these modes 
are contained in the RSC shell, which in all the cases have 10 layers. Also note in Fig. 
[TOT b) that the dipole is excited along the direction of the impinging wave, the x— axis. 
This property could be used as a sonic device to determine the direction of external sound 
sources and even to determine its exact positions by using an additional RSC shell[I]. In 
regards with quadrupole modes, a feature to remark is related to the additional number 
of oscillations that present the mode with frequency 0.6135 in comparison with that 
having lower frequency (0.5936). This difference is related with the higher number of 
oscillations that the periodic functions with larger argument Kq{Ra — Rb) have. 

5. Summary 

We have studied the resonances existing in two-dimensional (2D) RSC shells, which 
are acoustic structures consisting of a circular cavity surrounded by a finite RSC, both 
being embedded in a fluid background. Two type of resonant modes have been clearly 
characterized. On the one hand, cavity like modes, which have frequencies within the 
acoustic bandgaps of the RSC and are characterized by their strong localization inside 
the cavity. On the other hand, Fabry-Perot like modes, which appear in the frequency 
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Figure 11. Pressure maps (amplitude) of the first two Fabry-Perot like modes existing 
in the RSC shell with quadrupole symmetry (g =2). Their frequencies are 0.5936 
(a) and 0.6135 (b), respectively. They are excited by a punctual sound source (not 
depicted) put inside the central cavity; at {x,y) — {Q,d). The white circles define the 
borders of the RSC shell. 




Figure 12. Pressure maps (amplitude) calculated for the scattering of the RSC shell 
by an external sound with a plane wavefront impinging from the left side of the x— axis. 
The impinging sound has the same frequency than Fabry-Perot modes characterized 
in Fig. [nj(a) 0.5936 and (b) 0.6135, respectively. 
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regions where the RSC have an acoustic passband. The later modes are equivalent to 
those previously characterized in sonic crystal slabs[8] and are confined in the shell; i.e., 
between the cavity and the embedded background. 

The main advantage of the RSC shells is the possibility of selecting the mode 
symmetry either inside the cavity as well as in the shell by using bandgap engineering 
of the built in RSC. This mode symmetry selection can be used in designing devices 
for beam forming of the radiation field. Moreover, since the Fabry-Perot resonances are 
selectively excited along the direction of the external sound source, the dipole modes 
could be used in acoustic devices capable to determine the exact position of external 
sound sources. Moreover, RSC shells are excellent candidates for designing devices 
with functionality similar to the cochlea by following a recent proposal by Bazan and 
coworkers[9j, who used circular cavities with such purpose. In this regards, RSC shells 
have the advantage of their large tailoring possibilities. 

Let us remark that results obtained for 2D shells can be extrapolated to their three- 
dimensional (3D) counterparts, in which the cavity defect is a spherical hole. The high 
Q values that these structures present, either in 2D as well in 3D, make them suitable 
to fabricate acoustic resonators with possible applications in sensing and metrology [2] . 

With a few modifications, the results described here can be extended in the realm of 
electromagnetic waves, where radial photonic crystals (RPC) have been also proposed[l]. 
In this context, photonic cavities based on RPC shells are expected to have modes with 
huge Q— factors, a property that can be used to fabricate low threshold microcavity 
lasers with selected symmetry. 
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Appendix A. Demonstration of \Ms\ =1 
Let us define the matrix Ms(r) as 




(A.l) 



The matrix M, can be write as 



Ms = M;\Rb)Ms{Ra), 



(A.2) 



and the determinant of Mg is 




(A.3) 
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Applying the equation ([3D it can be shown that the derivative of determinant of Ms{r) 
is zero, so that this quantity is not a function of r, that is, 

dr\M,(r)\ = 



dr 



K{r)rpr ^{r)dr(f)„ (r) - dr (j)„ {r)rp^ ^(r)(9,.0+(r) 







so that 



(A.4) 
(A.5) 
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